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f) ' In axial symmetry, there is a gauge for Einstein equations such that the total 

^^, meiss of the spacetime can be written as a conserved, positive definite, integral 

on the spacelike slices. This property is expected to play an important role in 
kjq' the global evolution. In this gauge the equations reduce to a coupled hyperbolic- 

elliptic system which is formally singular at the axis. Due to the rather peculiar 
properties of the system, the local in time existence has proved to resist analysis 
^^ ' by standard methods. To analyze the principal part of the equations, which may 

represent the main source of the difficulties, we study linear perturbation around 
the flat Minkowski solution in this gauge. In this article we solve this linearized 
system explicitly in terms of integral transformations in a remarkable simple form. 
This representation is well suited to obtain useful estimates to apply in the non- 
linear case. 
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1 Introduction. 



The study of the gravitational systems with axi-symmetry is particularly appealing for at least 
5^ , a pair of important reasons. Firstly, quite a large number of interesting physical phenonrena are 

included and can described inside this category. Most notably the Kerr family and all its derived 
physics and mathematics belongs to it. Secondly because certain of its mathematical formula- 
tions enjoy a interestiirg number of remarkable mathematical properties. In particular there is a 
gauge, called the maximal-isothermal gauge in the system reduced by the axisymmetric Killing 
field where rather inrportant properties, as the positivity of nrass, are explicitly manifest and, 
presumably, would became important in the mathenratical investigations of these axisymnretric 
systems (see [I], [1], [3] and references therein). 



To take advantages of all this one needs to show, naturally, that the reduced Einstein equa- 
tions in the maximal- isothermal gauge is a mathematically well posed theory. As it usually 
happens in coordinates systems adapted to axial symmetry, the equations are formally singular 
at the axis. It is the particular combination of such formally singular terms inside the whole 
system of equations what makes it difficult to treat. Until now no well posedness result has been 
given in the literature. 

The analysis of the linearized equations inherit similar difficulties and therefore has proved 
to be non-trivial (see the discussion in [3], where this system was analyzed numerically). In 
this article we solve precisely this linear problem. The remarkable algebraic properties that this 
system has, allows us to solve it by a combination of integral transformations suitable adapted 
to these equations. This solution is the perfect analog to the solution in terms of Fourier 
transform of a constant coefficient equation. The construction appears to be finely adapted to 
this particular kind of singular hyperbolic-elliptic systems. 

We expect that this representation will be useful in the future to obtain relevant estimates 
for the non-linear case. In this sense, we believe that the present result opens the possibility to 
analyze the full axially symmetric Einstein equations in the maximal-isothermal gauge. 

The plan of the article is the following. In section [Ol we present our main result. At the end 
of this section we present the strategy of the proof, which is split in three main parts discussed 
in sections 12.11 12.21 and 12.31 respectively. Finally, we include an appendix in which we collect 
some properties of Bessel functions used in this article. 

1.1 Statement of the main result. 

In the maximal-isothermal gauge, the linearized Einstein equations with respect to a flat back- 
ground (in the twist-free case) reduce to the following system of equations for the functions v 
and /3 
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(2) 

See [1] for the deduction and physical discussion of these equationso 

In these equations the coordinates are (t, p, z). The relevant domain is the half plane < p, 
—oo<z<oo, which is denoted by M^. The axis is given by p = and it defines the boundary 
of the domain M^. A dot denotes time derivative, A is the flat Laplacian in 2-dimensions 

Av = dlv + a>, (3) 

and d denotes partial derivative with respect to the spatial coordinates p and z. 

The solutions {v{t, p, z), P{t, p, z)) we seek for in this article are C°° functions of M x M^ 
(meaning there is a C°° extension of {v, (3) into a open neighborhood of R x M"^ in R x M^). For 



^ We have slightly changed the notation with respect to this reference where the function /3 was denoted by 
pP to point out that it is the p component of the shift vector. Since in this article we will not discuss the other 
components of the linear perturbation (which can all be written in terms of v and j3) to simplify the notation we 
drop the index p. 



linearized gravity it is necessary to impose the following boundary and asymptotic conditions 
(see 15). 

For /3 it is required 

1. /3(i, 0, z) = and for every fixed (t, z), l3{t, p, z) is an odd function of p, 

2. For every fixed time (t), /3 = 0{r^^), and d^j3 = 0(r^^^'^), where r = y^p^ + z^. 
For V it is required 

1. v(t, 0, z) ~ {dpv)(t, 0, z) = and for every fixed (t, z), v{t, p, z) is an even function of p. 



2. For every fixed time (t), v = 0{r ^) and d'^v = 0{' 



,-2-fe 



Equation ((T]) is a wave equation for v and so it is neecesary to prescribe as initial data, 
roughly speaking, the position and velocity at the initial time (t = 0) which we will denote as 
u|t=o and i;|t=o- 

Remark 1. From the series expansion argument presented in OV it follows that, for analytic 
solutions, the requirement that v is even (in p) and (3 odd (for all times), follows only from the 
condition P{t, 0, z) = and u(0, 0, z) = {dpv){0, 0, z) ~ 0. 

The strategy to solve the system ([IJ-Q is to use an appropriate integral transformation to 
the whole set of equations to obtain simpler ones in the transformed variables. This integral 
transform is a combination of a Fourier transformation in the z coordinate and a Hankel trans- 
form in the p coordinate (see the appendix |^ for a definition of the Hankel transform) . The 
explicit form of the integral transform and its inverse are given by 



w(t, p, z) 

w(t, k, A) 
|A| 



w{t,k,X){k\X\p)iji{k\X\p)e'^'''^'dkdX, 
w{t,p,z){k\X\p)^Ji{k\X\p)e-^"^''dpdz, 



(4) 
(5) 



where Ji is the Bessel function of the first kind of order one. The ranges for the variables are 
< k < oo and — oo < A < oo. They define the same domain of integration as the variables {p, z) 
and hence we denote it by the same symbol K'^. By analogy to standard physical terminology 
we will call the space comprised by (p, z) the physical space, while the one comprised by (fc|A|, A) 
will be called the momentum space. 

Theorem 1. (Representation of solutions.) LetF{k,X) andG{k,X) be two arbitrary smooth 
functions of compact support in R1 . Define w and 7 by 



'w(t, k, A) 



fc2 



(f + fc2)i 7fc 



i^(/c,A)cos At(f + fc^)5 + 



and 



7(i,fc,A) = 2w{t,k,X) + 



2fct 

(f + P)f A 



G{k,X) 
A(f + fc2)- 

'(1 + fc')^ 

fc2 
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(Ai(f + P)^ 



w{t, k, A) dk. 



dk, (6) 



(7) 



For w and 7 define w and 7 by the integral transformation Q- Then, the functions v = p^w 
and /? = p~2 7 define a solution of equations H])-©, whose initial data f |t=o o,f^d i\t=o is given 
from F and G by 

v\t=o = ^P I , ^\. { r Ffk, A) d~k] (fc|A|p)^ Ji(fc|A|p)e2-^^ dkdX, (8) 

JmI (1 + k-^)2 \Jk J 

v\t=o = Vpf , ^' ,1 f /°°G(fc,A)dfcVfc|Ap|)^Ji(fc|A|p)e^"^-dfcdA. (9) 

Remark 2. /n this representation of solutions the boundary conditions required for j3 and v 
at the axis (items 1 above) are automatically satisfied. The asymptotic conditions required at 
infinity (items 2 above) are more difficult to establish, they require extensive technical analysis, 
and won't be included in the present article whose purpose is to introduce the generating formulas 
IQj and (^ in the most direct, simple and comprehensive fashion. In this sense we have preferred 
to use arbitrary functions F and G of compact support in the space (fc, A). This family of 
generating functions is at the same time simple and rich, and allows us to avoid technical lengthy 
developments in the proof of Theorem Ql 

The proof of theorem [T] is divided in three steps. In the first one, we use the scale invariance 
of the system to introduce new rescaled variables. The structure of the equations simplify in 
these new variables. This is done in section 12.11 The second step (section 12. 2p is to use the 
standard Fourier transform in the z coordinate. That essentially eliminates the z dependence 
of the equations. Finally, in section [^751 we analyze the p dependence of the equations using the 
Hankel transform. This is the most important part of the article. The proof of Theorem [T] is 
given thereafter. 

2 The equations in momentum space. 

2.1 Scaling symmetry 

Equations ((IJ-Q enjoy scaling symmetry (see [1]). This symmetry will play a fundamental role 
in the analysis of these equations. Let us describe this property. For a given solution v(t,p,z) 
and P(t, p, z) of equations ([I])-© we define the rescaled functions as 

vS,p,z)^v(i,PA), l3^^^/3(i,P/-), (10) 

\s s s J s \s s sj 

where 

t p z , ^ 

t=-, p^-, z = -. (11) 

s s s 

and s is a positive real number. Then, Vs and /3s define also a solution of equations ([T])-([2) 

written in terms of the rescaled coordinates (i, p, z). 

Note that v and /3 have different scaling. This difference manifests also in the power expansion 
series of these functions. 

To analyze the equations it is very convenient to introduce rescaled functions in such a way 
that both unknowns have the same scaling and the same behavior near the axis. Let us consider 



w and 7 defined by 

w = ^, 7 = Vp/3- (12) 

In terms of these variables, equations (dl)-© are given by 

zi) = Az.-^J + Vpa,(7p-^/^), (13) 



p A p^ \ 4 p- 



A7-^ + ,i=2 A.-^^). (14) 



The functions (^,7) are scale invariant in the following sense. For a given solution (^,7) of 
equations ([T^ - p^ the rescaled functions 

ftp z\ - ftp z\ 

w,s(t,p,z) =w [-,-,-] , 7s(i,p,z)= 7 -,-,-, (15) 

\s s sj \s s sj 

define also a solution of these equations in terms of the rescaled coordinates (fTTj) . Also, for a given 
solution, we expect w and 7 to have the same behavior at the axis, namely w = -f ~ 0{p^''^). 



Remark 3. The important part, of the rescaling (|12p is the relative power of p between v and (3 
which compensates the different scale behavior. We could consider scalings of the form v — p"w, 
P = p"~^7 for any arbitrary real number a. The specific power chosen in (J12p is motivated by 
the Hankel transform (see equations (J35p - (j36p in section\K^ . This choice makes the appropriate 
integral transformation symmetric with respect to its inverse. Another possible choice is a = 2. 
This power has the advantage that both functions w and 7 are even in p. However, the related 
Hankel transform is not symmetric with respect to its inverse and the formula does not coincide 
with the definition used in the literature. Hence to apply standard results concerning the Hankel 
transform we need to translate them into this new definition. This makes the proofs laborious, 
but there is no essential difficulty. 



Piw)^Aw-- — , (16) 



The differential operators in the spatial coordinates involved in this equation are given by 

3 w 

4^ 

and 

Th) = Vpdp {ip-'^') ■ (17) 

The distinction in the notation between P and T (boldface for P) is to emphasize that the 
differential operator P acts in both coordinates p and z while T does only in the p coordinate. 
Later on, we will define the operator P as the p part of P (see equation ([Ml)). As we will see, 
all the important features of the equations are contained in the p dependence. 

A remarkable fact is that the operators P^ and ^17} are also the natural operators for the 
second equation P^ . Namely, both equations ([T^ - P^ are written in terms of P and T as 
follows 

zi) = P(u;)+T(7), (18) 

P(7)-T(7) = 2P(^). (19) 

Note that this symmetry between both equations is not evident in terms of the original variables 
V and /3. The fact that in the second equation (|19p appears precisely this combination of P and 



T will be crucial. Our proofs will not work if we insert different (constant) coefficients in front 
of P and T in these equations. 

Remark 4. The operator P is singular at the axis. However, this kind of singular behavior is 
essentially the same as the one of the Laplace operator in n-dimensions for axially- symmetric 
functions written in terms of cylindrical coordinates (see, for example, the introduction of JB/)- 
A standard trick to avoid this problem is precisely to work in a higher dimensional space in 
which the operator is regular. This can be done also in the case of the operator P. However, the 
operator T will not be regular in that higher dimensional space. It appears not to be possible to 
write equations (|18p - (|19p as regular equations on a single higher dimensional space. 



It is the presence of the operator T in these equations which makes them so peculiar. The 
operator T is, outside the axis, a first order operator but at the axis it is a second order 
operator (due to L'Hopital rule). This behavior indicates that we can not decompose ([T5|- 
(|f 9p as a principal part (containing only the second order operator P) plus some lower order 
terms (containing only the operator T) . This kind of decomposition is essential to construct an 
iteration scheme in which each of the equations is solved in alternative steps of the iteration. 
Outside the axis this iteration scheme can be constructed, but it appears not to be possible 
to include the axis (see the heuristic discussion in |3]). In fact, our analysis suggests that the 
system (|18p -p ^ should be viewed as a unity that can not be further decomposed. 

2.2 The Fourier transform in z 

In equations (J18p - (|19p the z dependence is clearly simpler than the p dependence. The equations 
are regular in z and the coefficients of the differential operators do not depend on z. Hence, in 
order to factor out the z dependence we can use the Fourier transform in this coordinate defined 
in the standard way by 

/oo 
/(z)e-2-^^dz, (20) 

-OO 
/•OO 

r\f)^f[z)^ f{X)e^"^^dX. (21) 



Taking the Fourier transform in z to equations p^ -p ^ . we obtain the following equations for 
the transformed functions w{t,p,X), "f{t,p,X) 

{i,^Piw)-~X'^w + T{j), (22) 

Pi^)-X^^~Ti^)^2iP{w)-X^w), (23) 

where we have defined P as the p part of the operator P, namely 

Pi^)=^"-^- (24) 

A prime denotes derivative with respect to p. 

We use the scaling symmetry to reduce these equations to the case A = 1 in the following 
way. Define rescaled variables 

i=tX, p = p\X\, (25) 



then the rescaled functions wi(t,p), 7i(i,p) (no A dependence) satisfy the equations 

zi)i = P(wi)-wi+r(7i), (26) 

F(7i)-7i-T(7i) = 2(P(*i)-wi). (27) 



In these equations the derivatives are taken with respect to the rescaled coordinates ([25|) . If we 
have a solution wi{t,p), 7i(t, p) of equations (P5)) - (P7| . the solution of the original equations 
(Elll-dSSi) is given by 

w{t,p,X)=wi{tX,p\X\), 7(f,p,A)=7i(iA,p|A|). (28) 

The set of reduced equations (P5)) - (P7|) constitute our main equations. They encode all the 
main difficulties of the original equations. They will be solved in the next section. 

2.3 The Hankel transform in p 

To simplify the notation, let us write equations (P5 ]) -((77 |) without the tilde and without the 
subscript 1, namely 

w^P{w)-w + T{j), (29) 

Ph)-l-T{^)=2{Piw)-w). (30) 

The strategy to solve these equations is to expand the solution in terms of eigenfunctions of the 
operator P. That is, as a first step we look for a solution j of the eigenvalue equation for P 

P(i) = -fc2j. (31) 

By the same scaling argument used in the previous section (with A replaced by k), it is enough 
to consider the case k = 1 

Pi)) = -j. (32) 

Set j ~ y/pJ, then equation p2D in terms of J is given by 

p^J" + pj' + ip^ -1)J = 0. (33) 

This is the Bessel equation (see equation (|67p in the appendix). The behavior at the axis fixes 
the solution J to be, up to a factor, the Bessel function of the first kind of order one, denoted 
by Ji(jo). Thus, j = c^/pJi we take c = vfc. After rescaling. we have that the eigenfunctions of 
((3T|) are given by 

] = ^pJi{kp). (34) 

The solutions of equations (P5)) - ([5ni) will be constructed as a linear superposition of the eigen- 
functions j. As in the case of the Fourier transform with respect to the plane waves e^'^*'^^, the 
superposition of the eigenfunctions (|34p lead to the following integral transforms 



POO 

^(/) = />)=/ f{p)Vk~pJi{kp)dp, (35) 

Jo 

^-\f) - f{p) = / f{k)^pJi{kp)dk. (36) 







These are Hankel transformations of first order (see equations ([75|) - ([7il) and also [7] for further 
properties of the Hankel transform). The orthogonality property of the Bessel function (see 
equation (j69p l ensure that Sj~^fT{f) = / for / in an appropriate functional space (see the 
Appendix). We will call the space of f{p) the "physical space" and the space of f{k) the 
"momentum space". 

The rescaling P^ has been tailored to obtain precisely this form of the Hankel transform. 
Other rescaling will produce integral transforms which are not symmetric with respect to their 
inverse. They will have different weights in p. 

Let / = ^{f) and g = ^{g), then we have the following Parseval-type of identity 

OO /"OO 

f{p)g{p)dp^ / f{k)g{k)dk. (37) 

"'0 

That is, we have the following two identical inner products 

/•OO /"OO 

<f,9>^ f{p)9{p)dp, <f,g>= f{k)g{k)dk, (38) 

Jo JQ 

defined naturally in the physical and in the momentum space respectively. 

The crucial property of the Hankel transformation Sj for our purposes is its behavior with 
respect to the differential operator P, namely 

^(P(w)) = -k'^Sj{v). (39) 

This property follows after integration by parts and use of equation pip . Note also that the 
operator P is self adjoint with respect to < , >, namely 

< w, P{v) >^< P{w),v > . (40) 

2.4 The Hankel transform acting on the differential operators 

To solve equations ([^^ - ([5(I| we will apply the Hankel transform ([55)) to obtain simpler equations 
in momentum space. By equation (j39p . the Hankel transform i^ acts naturally on the operator 
P. However, this is not the case for the operator T. This operator is the main source of the 
difficulties. The next lemma characterize the action of Sj on T. 

Lemma 1. Let 7(fc) = -^(7) be of compact support. Then, the following relation holds 

^(T(7))(fc) = -kiE{k), (41) 

where we have defined 

E{k) = I %!^dk. (42) 

which is a function of compact support. 

Note that the function E{k) satisfies the following equation 

dE{k) _ 7(fc) 
dk ki 

This relation will be useful later on. 



(43) 



Proof. From the definition we iiave tiiat 

1"°° / \ 

i5(T(7))(fc) = y V^9p(7p-3/2J(pfc)^j,(fcp)dp, (44) 

and 

7(p)= / i{k){p~k)'^Ji{kp)dk. (45) 

Jo 

Multiplying by p^^'^ and taking a p derivative in equation (|45p under the integral, we obtain 

dp (IP-^A - rm'k^d, (:hM] dk. (46) 



We use the relation ([75]) to compute the derivative with respect to p of Ji . We have 

Mkp)_\' ^ ^,d / Ji(x)\ ^ _fc^^^^.^^_ ^^^^ 

/9 / ax \ X J p 

where we have defined x — kp. Hence we obtain 

dJ,p-i)^-rmk-'-^^dk. (48) 

^ ' Jo P 

Inserting the expression (|^ for dp {^p^^^"^) into equation (jii)) we get 

/•OO /"OO 

MTil)){k) = - k'^ki^{k)dk J2{kp).h{kp)dp. (49) 

Jo JQ 

From equation (j70p we know that the integral 

J2{kp).h{kp)dp, (50) 



is equal to fc/fc^ if A; < fc and it is equal to zero if A; > fc. Then, the conclusion of the lemma 

follows. n 

The next step is to analyze equation pO|) . In this equation there is no time derivative. We 
want to solve this equation for an arbitrary given function w, which is not necessarily a solution 
of the other equation (j29p . In the next lemma, we construct such solution using the Hankel 
transform. 

Lemma 2. Let w he a given function, with w{k) ~ Sj{'w) of compact support. Then, the solution 
7 = 55(7) of compact support of the equation (|30p in momentum space is given by 

7(fc) - 2w{k) + / ^ \, > wCk) dk. (51) 

We also have that E[k), defined by (|42p . is given by 

E{k)= ,^ ,,,. / ^ J. ' wCk)dk. (52) 

(1 + /C2)2 Jk fc2 



Proof. We apply the transform fj to equation ([50)1 . Using ([55)1 we get 

- (1 + fc2)^(fc) - i5(r(7)) - -2(1 + e)w{k). (53) 

We can use Lemma [T] to express S){T{'y)) in terms of 7. However, it is convenient to express 
everything in terms of E{k) (defined by (|42])) instead of 7(fc). Using (|41]), (|42]) and (|43]) we 
obtain 

We multiply by the integrating factor (1 + fc^)^ to get 

A((i + fc2)i^)^_2ii±pizi(fc). (55) 

Integrating this equation and forgetting the integrating constant to have a solution of compact 
support, we obtain equation ([5^ . Equation (|5T|) follows directly using (1^5)) . D 

Finally, in the next lemma we solve the whole system ([^^ -([50 )) . 

Lemma 3. Let F{k) and G{k) be arbitrary functions of compact support. Then w = Sj^^{w) 
and 7 ~ 53^^(7) define a solution of the system of equations (J29|) - p0p . where 

w(t,k) = (F(k)cos((l + P)h) + G(k)sin((l + P)h)]dk, (56) 

(l + fc2)2 7fc ^ ^ 

j{t,k) = 2w{t,k)+ ^ T., ^ w{t,k)dk. (57) 

Proof. First, note the equation ([57)1 is the solution of equation ([50]) given by lemma [5] if we 
consider w(t, fc) as a given function. The only part we have to prove is that equation (|56p is also 
a solution of (|29|) in which 7 is given by (|57)l . 

We apply the Hankel transform to equation (P^)) and we obtain the following equation in 
momentum space 

ib{t, k) = -(1 + fc2)u,(t, fc) + io(T(7)). (58) 

Using the expression for S){T{'-^)) obtained in lemma [2] we obtain 

w{t, k) = -(1 + e)w{t, k) - / ^ \, '' wit, k) dk. (59) 

(1 + k^)2 Jk k2 

This equation involves only li, and hence we have reduced the system ([^^ -([50 ]) to only one 
equation for one unknown. But equation (j59p is not a differential equation, it is an integro- 
differential equation. We can get a simpler expression if we define a(t, k) by 

a{t, k) = ^ 3 '' w{t, k). (60) 

fc2 

In terms of a{t, fc), equation (|59p is written as 

/>oo 

a(i,fc) = -(l + fc2)a(t, fc)-2 / ka{t,k)dk. (61) 
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Although equation (pT|) looks certainly simpler than equation (|60p , the essential difficulty remains 
the same. However, a remarkable cancellation occurs if we take a derivative with respect to k 
allowing us to convert this equation into a pure differential equation. Namely, let us define 

b{t,k) = ^a{t,k). (62) 

Then, taking a derivative with respect to k of equation (|6ip we obtain the following remarkable 
simple differential equation for b(k) 

b{t,k)^-{l + k^)b(t,k). (63) 

The important fact is that the term —2ka{t, k) that appears when the k derivative is applied to 
—k^a{t, k) in equation ((CT|) cancels out by the derivative of the integral (the factor 2 in front of 
the integral is crucial). And hence in the final expression only appears 6(i, k) and not a(i, k). 
The solution of equation (j63p is given by 

6(t, k) = -F(fc) cos((l + k^)^t) - G{k) sin((l + k^)^t), (64) 

for arbitrary functions F{k) and G{k). We have written equation (|64p with a minus sign just 
for convenience. 

The function a{t,k) is calculated integrating (|62p . namely 

a(i,fc) = - / b{t,k)dk. (65) 



There is a constant of integration that we have set to zero in (j65p . otherwise the function w 
defined by (pO|) will not be of compact support. 

Using dni), dMI) and ([501), expression ^ follows. D 

3 Proof of Theorem [H 

Proof. The theorem is a straightforward consequence of lemma [3l and the scaling invariance of 
the equations. Namely, consider the solution of equations ([^^ - ([50]) founded in lemma [3] To be 
consistent with the notation used in section [521 this solution should be denoted by 'Wi{t,p) and 
ji{t,p). Using the scaling ([25)) . out oi wi{t, p),'yi{t,p) we construct 

w{t,p,X)^wiiXt,\X\p), 7(<,p,A)=7i(Ai,|A|p). (66) 

These are solutions of (P^ - (P^ . We apply the inverse Fourier transform ([51]) to ■w{t,p,X) and 
j{t,p,X) to obtain the desired result. In © we have redefined the function G{k,X) to make 
simpler the connection with the initial data. D 
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A Bessel functions and Hankel transform 

We collect in this appendix properties on Bessel functions and Hankel transform that we use in 
this article. A general reference is the classical book and also the book [7| for the Hankel 
transform. 

We will consider Bessel functions of the first kind, denoted by J^ (x) . We will restrict ourselves 
to the case where i/ is a positive integer. 

Bessel functions are solutions of the Bessel equation 

''■'"■i^+(l-^)^.^0. (.7) 



dx^ X dx \ X _ 
For i^ > 0, a series expansion of these function is given by 



where F denotes the standard Gamma function. Since in our case i^ is a positive integer we have 
r{j + V + 1) = {j + v)\. From ([55)1 we deduce that J,y{x) is an even function of x for v even and 
an odd function of x for i/ odd. 

We have the following orthogonality property for Bessel functions 

1 - 
pJJkp)JJkp)dp = -5{k - k). (69) 

k 

The following integral is used in the proof of lemma [1] (see [S] , p. 406) 

J^{at)J^_i{ht)dt = { 1/(26), (70) 

.0, 

where the first value corresponds to & < a, the second to 6 = a and the third to 6 > a. 
We make use also of the following relation ((5] p. 45) 

X — -^ vJ^{x) ^ -xJ^+i{x). (71) 

ax 

For the case v — \ this relation can be written in the form 

d ( Ji{x)\ _ J2{x) 



,,, . , . m 
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The Hankel transform of order v of the function / is defined as 



^.(/) = / f{pWkpUkp)dp, (73) 

and the inverse is given by 

^Z\F) = / F{k)^/k^Mkp)dk. (74) 

The formula (|69)) guarantee that i5~^(io^(/)) — f with / in an appropriate functional space. 
More precisely the Hankel transform (of order one) is an automorphism on the space Hi oi 
C°° functions / : M+ -^ M provided with the family of seminorms 7m, fc which for each pair of 
non-negative integers m and k are defined as 

7™,fc(/)= sup \p"\p-'d,f[p-''^f{p)]\. (75) 

0</9<OO 

A function / in Hi has an expansion of the form 

f{p) = pi{ao + a2P^ + ... + a2ip^^ + R2i{p)), (76) 

and is of fast decay, namely, decays faster than any power of 1/p as p — > oo (Lemma 5.2.1 in pg 
130 of CD). 
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